Height representation, critical exponents, and ergodicity in 
the four-state triangular Potts antiferromagnet 
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We study the four-state antiferromagnetic Potts model on the triangular lattice. We show that 
the model has six types of defects which diffuse and annihilate according to certain conservation 
laws consistent with their having a vector-valued topological charge. Using the properties of these 
defects, we deduce a 2 -I- 2-dimensional height representation for the model and hence show that 
the model is equivalent to the three-state Potts antiferromagnet on the Kagome lattice and to 
bond-coloring models on the triangular and hexagonal lattices. We also calculate critical exponents 
for the ground state ensemble of the model. We find that the exponents governing the spin-spin 
correlation function and spin fluctuations violate the Fisher scaling law because of constraints on 
path length which increase the effective wavelength of the spin operator on the height lattice. 
We confirm our predictions by extensive Monte Carlo simulations of the model using the Wang- 
Swendsen-Kotecky cluster algorithm. Although this algorithm is not ergodic on lattices with toroidal 
boundary conditions, we prove that it is ergodic on lattices with free boundary conditions, or more 
generally on lattices possessing no non-contractible loops of infinite order. To guard against biases 
introduced by lack of ergodicity, we therefore perform our simulations on both the torus and the 
projective plane. 



I. INTRODUCTION 

There has been considerable interest over the last few 
years in the properties of classical spin systems possess- 
ing highly degenerate ground states. Many such models, 
including ice models [y|2|, the triangular Ising antiferro- 
magnet [pUl], and dimer models MM have been found 
to have ground state ensembles which display critical 
properties such as algebraically decaying spin-spin cor- 
relations and divergent fluctuations in the order parame- 
ter. It is now known that these properties are associated 
with the existence of interface or "solid-on-solid" repre- 
sentations for the models, in which sites can be assigned 
heights h which vary smoothly over the lattice and which 
can be mapped onto the states of the spins or other mi- 
croscopic variables in a simple way. If we assume that 
this height field behaves as a Gaussian surface, which is 
justified if the model is in its rough phase, then the crit- 
ical behavior follows in a straightforward fashion [0-^. 
The values of the critical exponents are related to the 
stiffness K of the surface and the wavelength of the ap- 
propriate operator on the height lattice. 

A few models have been studied for which h is vector 
rather than scalar and the ideas above generalize to this 
case also M-Oj. In this paper we look at one particu- 
lar model of this type, the four-state antiferromagnetic 
Potts model on the triangular lattice. This turns out to 
be an especially lucid example of a model with a vector 
height, being defined on a simple Bravais lattice and, as 
we will show, possessing a very straightforward height 
representation. 

It is clear that the four-state antiferromagnet does in- 
deed have a highly degenerate ground state, since the 



triangular lattice is three-colorable. By taking a three 
coloring and introducing a finite density of the fourth 
color we can see that the ground state must have an ex- 
tensive entropy. In this paper, we study the properties of 
the ground state ensemble by considering first the behav- 
ior of defects in the model at zero temperature. We show 
that there are six distinct types of non-trivial defects and 
from the conservation laws that govern their collisions we 
deduce that they have vector charges ^tt apart. We use 
this observation to derive a Burgers vector for the model 
and hence show that when no defects are present the sys- 
tem has a two-dimensional height representation. The 
defects then correspond to screw dislocations on a 2 -f 2 
dimensional lattice and we predict that pairs of them will 
be attracted or repelled with an entropic Coulomb force 
proportional to the dot product of their charges. 

We use our height representation to deduce a num- 
ber of facts about the four-state model. First we show 
that it is equivalent to the three-state Potts antiferro- 
magnet on the bonds of a hexagonal lattice. This equiv- 
alence has also been derived using a different approach 
by Baxter uM, but the derivation given here is nonethe- 
less instructive because it respects the symmetries of the 
system under permutation of states in a way that Bax- 
ter's does not. By a simple geometrical construction we 
show further that the model is equivalent to the g = 3 
antiferromagnet on the Kagome lattice, a model which 
has previously been studied by Huse and Rutenberg [10 
And employing results due to Kondev and Henley [15 
we show that our model is equivalent to the fully-packed 
loop model on the honeycomb lattice with a loop fugac- 
ity of 2. We generalize these equivalences to several other 
cases, including one related to loop models on the square 
lattice. 



The existence of a height representation also impUes, as 
mentioned above, that the ground state ensemble is crit- 
ical and we have verified this by Monte Carlo simulation. 
Simulation of this model is not trivial, since no single- 
spin-flip algorithm is crgodic and the best-known cluster 
algorithm, that of Wang, Swendsen and Kotecky ||l6| , |l7| , 
is believed not to be ergodic under toroidal boundary 
conditions [Q and has not been proved to be ergodic 
for any other case. (For other models, however, it is 
known to be ergodic, particularly for models defined on 
bipartite lattices p^ , p^ , pO| |.) Here we make use of our 
height representation to prove for the first time that the 
algorithm is in fact ergodic for the q = A Potts model 
on any lattice with triangular plaquets satisfying certain 
topological conditions. This includes lattices with free 
boundary conditions or with the topology of a sphere or 
a projective plane. To take advantage of this result we 
have performed our simulations on the projective plane. 
The idea of changing the lattice's large-scale topology in 
this way to make a sampling algorithm ergodic appears 
to be new. 

The outline of the paper is as follows. In Section Owe 
introduce the model and study the types of defects which 
occur in it and their interactions. In Section III we de- 



rive the height representation of the model and thereby 
demonstrate the model's equivalence to various others. 
In Section f^ we show how the scaling exponents charac- 
terizing the large-r behavior of correlation functions can 
be calculated from the height representation and in Sec- 
tion |v| we demonstrate the existence of entropic Coulomb 
forces between defects. In Section VI we describe the 



Monte Ca rlo a lgorithm we use to simulate the model and 
in Section |VII| g ive the results of our simulations. Finally, 
in Section [VI 1 1| we give our conclusions. 



II. DEFECTS IN THE FOUR-STATE 
TRIANGULAR POTTS ANTIFERROMAGNET 

The g-state Potts model is a generalization of the Ising 
model in which a lattice is populated with spins Si , one on 
each vertex i, which can take integer values Si — 1 . . .q. 
The spin states are also sometimes referred to as "colors," 
and we will occasionally make use of this metaphor. The 
energy of a configuration is defined to be proportional to 
the number of pairs of adjacent sites with the same state 






(1) 



In this paper we consider the Potts model with q — A 
on the triangular lattice in two dimensions and with 
J < which makes the model antiferromagnetic so that 
similarly-colored pairs of adjacent sites are energetically 
unfavorable. We refer to such pairs as "defects." This 
model can also be thought of as a discretization of the 
classical Heisenberg model in which each site has a three- 
dimensional unit vector spin S^: 
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FIG. 1. Defects with the same state on both sides of the 
frustrated bond we call "false," because they can be annealed 
away without interacting with any other defects. 
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This is equivalent to (M up to a rescaling and an additive 
constant if the Sj are restricted to the corners of a tetra- 
hedron, since then the dot product of two spins depends 
only on whether they arc the same or different. 

The ground state entropy per site for this model can 
be calculated analytically |2^ or closely bounded with se- 
ries approximations [^; its exact value is 3r(i) /47r^ ~ 
1.460998. 

Consider the behavior of the model under a single-spin- 
flip dynamics at zero temperature. Such a dynamics al- 
lows existing defects to diffuse and interact in a variety 
of ways, but creates no new defects. We now show that 
the defects in the model fall into a number of different 
classes with well-defined properties. The simplest case is 

b 
that of a defect of the form a a in which the sites to 

b 
either side of the defect pair both have the same state 
(here denoted b) . As Figure |l| shows, these defects retain 
this same form when they move. Moreover, it is possible 
for a defect of this type to disappear entirely if it en- 
counters a neighborhood of the right configuration. This 
is shown in the rightmost portion of the figure, where 
the center site in the hexagon can become a 4, resulting 
in a defect-free configuration. Since these defects can be 
annealed away with only local moves, and do not require 
interaction with any other defects, we call them "false" ; 
their density can be expected to fall off exponentially fast 
in a quench to T = and so can be ignored where the 
long-time relaxation of the system is concerned. 

b 

This leaves us with defects a a whose neighboring 

c 
sites are in different states b ^ c. As Figure shows, 
these defects are persistent and an isolated defect cannot 
be annealed away by local moves alone. The rest of our 
analysis in this section will concentrate on these "true" 
defects. 

The true defects possess two properties which are con- 
served both during diffusion and in interactions with 
other defects. The first depends on the states a, b and 
c which make up the defect and its immediate neighbor- 
hood. Since a, b and c are all different for an isolated true 
defect, they divide the four spin states in the model into 
two pairs, where a belongs to one pair, and the sites on 
either side belong to the other. In the case of the defect 



23 
1 1 4 

3 2 



3 
-44 

2 

24 

1 1 2 

3 4 



23 

1 1 2 

3 4 

.1 3 

3 4 



42 
34 

24 

1 1 3 

3 4 



24 

1 1 3 

32 

. 24 
1 2 



24 
43 



23 



2 1 



2 1 



X+X-0 1 1 1-0 1 1 3-0 1 1 3 




FIG. 2. As a defect diffuses, it stays witliin one of six equiv- 
alence classes, as defined in the text. Here we show the inter- 
action of a defect of type X with all possible neighborhoods 
(excluding symmetry equivalents). 
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FIG. 3. The defect classes X, Y and Z all correspond to 
counter-clockwise movement on the outside of one of the faces 
of the pictured tetrahedron. The defects X, Y and Z are their 
mirror images and correspond to clockwise movement. 



1 1 , for example, the defect sites belong to the pair 

3 
{1,4} and the adjoining sites to the pair {2,3}. There 
are three distinct ways of dividing up the states in this 
fashion. 

The other conserved property of a true defect is a 
handedness defined as follows. The states a, b and c, 
in that order, describe a path which is either clockwise 
or counter-clockwise on the outside of one of the faces of 
a tetrahedron whose vertices are labelled with the four 
spin states as shown in Figure 0. In order that this prop- 
erty be correctly conserved we must in addition stipulate 
that it remains the same under 120° rotations of a defect, 
but changes sign under 60° or 180° ones. Thus for ex- 

2 
ample the defect 1 1 considered above has a counter- 

3 
clockwise (or positive) handedness, while an inversion or 
60° rotation gives us a clockwise (or negative) handed- 
ness. 

Using these two conserved properties, we can divide 
the 72 possible defects into 6 equivalence classes which 
we label X, Y , Z and X, Y , Z. Here the letters denote 
the pairing of states and the bars (or absence of them) 
denote the handedness. Representative members of the 
classes are shown in Figure 0. 
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FIG. 4. Examples of collisions between defects. Note that 
in collisions such as X -I- X —> 0, a false defect is construed as 
no defect at all. 




FIG. 5. The vector charges assigned to the defects. The 
vectors a, (3 and 7 are defined in the text. 



Considering again a single-spin-flip dynamics, we show 
in Figure a selection of possible collisions between var- 
ious types of defects. Although more complex collisions 
than these can occur, we always find (and it is proved 
below) that X +li ^ 0, X + Y -> Z, and similarly 
for cyclic permutations. If we wish to assign a charge 
X to each particle such that x(^) + xO^) + x{Z) — 0, 
x{X) = — x(X), and so on, we can do this with six vec- 
tors ^TT apart as in Figure g|. We adopt the convention 
that Ixl = 1, X(^) = (1,0), and x{Y) = (-1/2,^3/2). 
The proof that charge is locally conserved goes as follows. 

It turns out that the total charge within an area 
can be expressed quite simply as a sum over the pla- 
quets making up that area. We start by writing the 

a 
charge inside a diamond he as a sum of functions 

d 
of the upward- and downward-pointing triangular pla- 
quets, /(a, 6, c) +g{d, c, b). Since a 180° rotation reverses 
the charges, we have g = — /, and since a 120° rotation 
keeps them the same, / is symmetric under cyclic permu- 
tations, /(a,6, c) = f{b,c,a) = f{c,a,b). Furthermore, 
since a defect- free plaquet has zero charge, f{a,b,c) — 
if a, 6,c are all different, and since rotating the tetrahe- 
dron of Figure |3| around the corner corresponding to spin 
state a rotates the charge plane of Figure |^ around the 
origin, f{a, a, a) = for any a. 

This just leaves the case where exactly two of a, b and 
c are equal. Solving the equations /(2, 1, 1) — /(3, 1, 1) = 
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FIG. 6. Left: the direction of the vector t along 
the edges of the triangular lattice. Right: the sum 
j> ds ■ B — a ~ P = 1 around a diamond enclosing an X 
defect and a/2 + /3/2 + 7/2 = around a defect-free triangle. 
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FIG. 7. Left: the mapping from heights to states. Right: 
integrating around a defect, in this case an X, gives a Burgers 
vector of A/i = %/3 (2a - 2/3) = 2^/3 x- 



X{X), /(3, 1, 1) - /(4, 1, 1) = x{Y), and so on, gives 



fia,a,b) = 



if {a, fe} 
if {a, b} 
if {a, b} 



{1,2} or {3,4} 
{1,3} or {2,4} 
{1,4} or {2,3} 



where 



a 



f3 



7 



2' 2^/3^' 

2'2V3/^ 
M 



(°'-;7i 



V3^' 



(3) 

(4) 
(5) 
(6) 



as shown in Figure pi 

If the charge is conserved in all collisions as we have 
claimed, we should be able to write it as an integral of 
some quantity around the perimeter of the area we are 
interested in. It turns out that this is indeed possible if 
we define a tensor B on each edge of the lattice equal 
to the outer product B = t (g) E of two vectors, t and 
E. The first of these points along the edge and gives 
it a direction as in Figure g. This ensures that B has 
the necessary change of sign under 60° rotations. The 
second is a vector in the charge plane which depends 
symmetrically on the states a and b at the two ends of 
the edge thus: 



E 




if {a, 6} -{1,2} or {3, 4} 
if {a, 5} = {1,3} or {2, 4} 
if {a, 5} -{1,4} or {2,3} 
a a — b. 



(7) 



Then the charge inside a finite region of the lattice is an 
integral around a counter-clockwise perimeter 



(b ds-B = ^(As-t)E. 



(8) 



In Figure |6| we show this for a diamond around a sin- 
gle defect and for a defect-free triangle. Since all larger 



regions can be formed by gluing together diamonds and 
triangles like these, and since the integral cancels along 
the shared edges within the region because of the sign 
change imparted on i? by t, it follows that the charge 
within any region is correctly given by Equation (0). This 
proves our contention that the charge is conserved in all 
defect collisions, since the value of the integral around 
a line completely enclosing any such collision does not 
change when the collision takes place. 



III. HEIGHT REPRESENTATION AND 
EQUIVALENCE TO OTHER MODELS 

The integral in Equation (0) defines a Burgers vector 
for a defect in our model. Around any defect free re- 
gion, the Burgers vector is zero, and hence on a lattice 
possessing no defects the integral of B between any two 
points is path independent. This allows us to define a 
height representation for the model in which the height 
h at any site is specified uniquely as the integral of ds • _B 
from a single reference site of known height to the site 
of interest. The heights are thus, like the Burgers vector 
itself, two-dimensional vectors living on a triangular lat- 
tice with lattice vectors a, f3 and 7, multiplied by ±V3 
in order to make the lattice constant 1. In fact, it is 
straightforward to show that there is a unique mapping 
of heights onto spin states, which is a four-coloring of 
the height lattice as shown on the left-hand side of Fig- 
ure M. The particular permutation of the colors in this 
figure depends on the definition oi a, (3 and 7 given in 
Equations (m through (H) and on the choice of reference 
site. 

Once we have the height representation for the model, 
there are a number of results which follow. In this sec- 
tion, we use it to demonstrate the equivalence of the 
ground state ensemble to a number of other models, some 
of which have been studied previously. 

First, imagine coloring the edges of the triangular lat- 
tice in a defect-free configuration of the model with three 
colors a, /3 and 7 according to the height difference along 




FIG. 8. The relationship between the four-state antiferro- 
magnet on the triangular lattice (dashed) and the three-state 
antiferromagnet on the vertices of the Kagome lattice (bold) 
or the edges of the honeycomb lattice (dotted). 

them, as in Figure |[ If we define two edges as neighbor- 
ing when they bound the same triangle, then neighboring 
edges must have different colors since otherwise two of the 
vertices of the triangle would have the same spin state. 
Thus the model is equivalent to a three-coloring of the 
bonds of the triangular lattice in which no two adjacent 
bonds are the same color. 

The reverse mapping is also possible. Since a + P + j = 
0, the change in height sums to zero around any triangle, 
and therefore around any closed curve, so the height and 
therefore the state of every site is well-defined once we 
have chosen the spin state of one reference site. As there 
are four choices for this reference state, every configu- 
ration of this bond-coloring model corresponds to four 
ground states of the four-state triangular Potts model. 

Since the edges of a lattice are in one-to-one correspon- 
dence with the edges of its dual lattice, we can also think 
of the model as a three-coloring of the edges of the honey- 
comb lattice, where two edges are neighbors if they share 
a vertex (see Figure |8|) . A simple extension of these map- 
pings is to put a vertex at the midpoint of each edge on 
the triangular lattice (or the honeycomb lattice) and con- 
nect those vertices which fall on neighboring edges. The 
result, as shown in Figure ||, is a Kagome lattice. (In 
general, this construction is called the "medial graph" 
|P3| .) Thus the four-state triangular Potts antiferromag- 
net and the three-state one on the Kagome lattice also 
have equivalent ground-state ensembles. 

A number of these results have appeared previously 
in one form or another. Huse and Rutenberg pO| found 
a two-dimensional height representation for the q — 3 
antiferromagnet on the Kagome lattice, which is equiva- 
lent to ours once the equivalence between models demon- 
strated above is taken into account. Baxter m^ (see 
also Ref. [^ ) demonstrated the equivalence of the four- 



state antiferromagnet and the bond-coloring model on 
the hexagonal lattice using an approach somewhat dif- 
ferent from ours. He defined a cyclic ordering 4 — > 3 ^ 

2 — » 1 — > 4, drawing arrows from higher states to lower 
ones (modulo 4) and leaving edges between states 1 and 

3 or 2 and 4 blank. However, since the four-state an- 
tiferromagnet is invariant under all permutations of the 
four states, not just cyclic ones, we feel that the mapping 
given here better respects the symmetries of the system 
under such permutations. 

Kondev and Henley psf showed that the bond-coloring 
model on the honeycomb lattice is also equivalent to a 
fully packed loop (FPL) model on the honeycomb lat- 
tice, where loops are defined as sets of edges alternating 
between two of the three colors. The loops are then con- 
tours of the component of the height perpendicular to 
the direction corresponding to the third color [g5| . Each 
loop can have its colors exchanged without affecting the 
surrounding configuration. Such loops arc said to have a 
fugacity n = 2. (When n — 1, the FPL model is equiv- 
alent to the triangular Ising antiferromagnet |26|.) An 
exact solution for the ground state entropy of the FPL 
model on the honeycomb lattice for general n has been 
given by Batchelor, Suzuki and Yung p^ using a Bethe 
ansatz, and differs from Baxter's solution for the entropy 
of the q — A triangular Potts antiferromagnet |gl| by 
exactly log 4, as we would expect given the equivalence 
demonstrated above. 

The equivalence between four-state antiferromagnets, 
three-state bond-coloring models, and fully-packed loop 
models applies on other lattices as well. If a lattice has 
triangular plaquets and vertices with even coordination 
number, we can define an orientation on the bonds with 
vectors t which go either only clockwise or only counter- 
clockwise around each plaquet. Then we can define B — 
t(g)E as before, and any path around a defect-free plaquet 
will have Ah = zL{a + P + ^) = 0, so that a consistent 
definition of h exists. Furthermore, if we color the bonds 
with colors a, f3 and 7 then each vertex of the dual lattice 
has one bond of each color. Bonds of any two colors 
comprise a set of fully packed loops with fugacity 2, and 
these are contours of the component of h perpendicular 
to the remaining color. 

As an example, the lattice shown on the left-hand side 
of Figure has vertices of coordination number 4 and 8, 
and the orientation of the bonds can be defined as shown. 
The q = A Potts antiferromagnet on this lattice is equiva- 
lent to a bond-coloring or fully-packed loop model on its 
dual lattice, the truncated square lattice, whose plaquets 
are squares and octagons. If we take the trace over the 
sites with four neighbors, we are left with a model on 
the square lattice where plaquets of the form ° , are 



prohibited, while those of the form 



and 



have 



fugacities 1 and 2 respectively. This model was studied 
by Nienhuis ||2^ and is also equivalent to a model on the 
square lattice where loops can collide but not cross. A 
similar model where the latter two types of plaquet have 
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FIG. 9. Left: another lattice on which the four-state anti- 
ferromagnet has a height representation (hnes with arrows), 
and the corresponding dual lattice (dotted). Right: the 
fully-packed loops of a and f3 bonds corresponding to the 
spin configuration shown. This model is also equivalent to a 
loop model on the square lattice where loops can collide but 
not cross. 



equal fugacity was studied by Burton and Henley [[l3[ , 
who found a five-dimensional height representation for 
it. 

We close this section by defining another model equiv- 
alent to the (7 = 4 antiferromagnet. Let us define a chi- 
rality on each plaquet of the triangular lattice according 
to whether the states at its three corners are oriented 
like an interior or exterior face of the tetrahedron in Fig- 
ure 1^, or equivalently, whether the colors on its edges in 
the three-coloring model above cycle clockwise or anti- 
clockwise. These plaquets correspond to the vertices of 
the honeycomb lattice, and it is not hard to see that 0, 
3, or 6 of the vertices of each hexagon must have positive 
chirality. There are 12 configurations of the four-state 
triangular antiferromagnet for each state of this model, 
one for each choice of the states of two adjacent reference 
sites. 



IV. FREE ENERGY AND CALCULATION OF 
SCALING EXPONENTS 

Consider the restricted entropy of the four-state trian- 
gular model for a particular value of some (unspecified) 
coarse-graining of the height field h. It is not hard to 
convince oneself that this entropy is lowest when |V/i| is 
large, and highest for configurations that are macroscop- 
ically flat. For instance, a three-coloring of the lattice is 
flat since any two sites with the same spin state have the 
same height, and the set of configurations in the vicinity 
of such a three-coloring contributes a large entropy to 
the ground-state ensemble since every site has a choice 
of two colors. On the other hand, a four-coloring whose 
height increases linearly across the lattice corresponds to 
only one microstate, since no site has any choices at all. 
Building on considerations such as these we can derive 
expressions for the scaling exponents of the model. Our 
presentation follows that of Burton and Henley E3| . 



If the model is in its rough phase, the arguments of the 
previous paragraph suggest that it has an effective free 
energy of the form 



G = -K"^ 
2 



S/h^\/hxdxdy, 



(9) 



where i^"^ is a stiffness tensor. However, since the model 
is invariant under permutations of the four spin states, 
and since these permutations are equivalent to rotations 
of the height lattice, K must be a scalar and we have 



2 



(|V/ii|2-K|V/i2p) 



dxdy, 



(10) 



where hi and ft,2 are components of the height field along 
any two perpendicular directions in the height space. (In 
our calculations we have taken hi and /i2 along the direc- 
tions indicated in Figure M.) In frequency space this free 
energy decouples into a sum over independent Gaussians, 
giving spatial correlations between the heights at points 
a distance r — |r2 — ri| apart of 



{\h{r,)~h{ri)\')^-Liogr + C 



(11) 



for large r, where C is a constant. 

Quantities such as spin and local magnetization are 
periodic functions f{h) of the height and hence can be 
Fourier expanded in the height. We can calculate the 
spatial correlations in any one such Fourier component 



f9 



cx e 



ig-h 



having frequency g in height space, using 



the fact that h's Fourier components are Gaussianly dis- 
tributed. This gives us 

(/g(ri)/,(r2))«(e'^-l''('-^)-''('-^)l) 

-p-^/{l'>(r2)-?i(ri)|^> 



-(d-2+ri) 



(12) 



for large r where rj is the anomalous dimension of the 
correlation function and d is the dimensionality of the 
lattice. Given that d = 2 in the present case and making 
use of Equation (O) we then find that 



V = 



9 



2TrK 



(13) 



If a quantity has several non-zero Fourier compo- 
nents, then the one with the smallest rj — i.e., longest 
wavelength — will dominate for large r. 

If K is not a scalar, these equations generalize to 



{{h',-h,){h'x~-hx)) 



k: 



■ log r -I- Ck> 



and 



27r 



(14) 



(15) 



For instance, transforming a scalar K to triangular co- 
ordinates where basis vectors are jjtt apart rather than 
orthogonal gives a matrix K' 



K' ^ K 



1 

-1/2 



-1/2' 



and {K' 



J_ /4/3 2/3 
if U/a 4/3 



(16) 



Using the equivalence of the bond-coloring model to 
the fully-packed loop model on the honeycomb lattice, 
and relating the vortex-antivortex correlation function 
to the probability that two sites lie on the same loop, 
Kondev and Henley [^ have shown that these models 
are exactly at their roughening transition and have a stiff- 
ness oi K = irTT. Since our four-state model is equivalent 
to these models, it has the same value of K. We will use 
these results below to calculate the scaling exponents of 
various quantities for comparison with Monte Carlo ex- 
periments. 



V. FORCES BETWEEN DEFECTS 

Since the energy of a pair of defects is 2 regardless 
of how far apart they are, there is no energy gradient to 
drive a force between them. However, there is an entropic 
force, driven by the fact that the presence of a free defect 
reduces the entropy within an area of radius r by an 
amount proportional to logr. 

In a model with a one-dimensional height representa- 
tion, a defect with Burgers vector b has an average field 
around it 
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b 
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(17) 



giving a force between two defects with Burgers vectors 
b and b' 



F 



Kbb' 



(18) 



When coupled with a mobility F, this gives an average 
velocity to the defects of 



TF 
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(19) 



Such forces have been measured numerically for the 
three-state Potts model on the square lattice [Q. The 
generalization to a higher-dimensional height representa- 
tion is straightforward. Since the free energy in Equa- 
tion (|lO|) is a sum of independent terms in hi and /12, the 
force between two defects with Burgers vectors b and b' 
will be proportional to bib'i -\- 62^2 = b • b'. Since the 
Burgers vector for a defect is equal to 2-\/3 times the 
charge x on that defect (see Section HI), this gives us a 
force of 



F 



12Kx-x' 



(20) 



In other words, an X and a Y will be attracted to each 
other, but only half as strongly as an X and an X, and 
an X and a Y will be repelled half as strongly as an X 
and another X. 



VI. MONTE CARLO SIMULATION 

In order to simulate correctly the properties of the 
ground-state ensemble of a system, it suffices to find a 
set of update moves which take us from one ground state 
to another without introducing any defects, such that 
every ground state can be reached from every other in a 
finite number of moves on a finite lattice. An algorithm 
based on such a set of moves is said to be ergodic, and 
it can be shown that any ergodic algorithm will sample 
all ground states with equal frequency over the course of 
a long simulation (see Ref. [Q for example). Unfortu- 
nately, it turns out to be quite difficult to find a suitable 
ergodic set of moves for the four-state Potts model con- 
sidered here. To begin with, it is clear that no single- 
spin-flip dynamics can be ergodic because finite defect- 
free regions can be pinned under such a dynamics, and 

remain pinned no matter what happens outside them. 

1 2 

For instance, in the hexagon 343 every spin has at 

2 1 

least one neighbor of each of the other states, so no spin 
can change state. Since there is a positive density of 
such clusters in a random ground state on a large lattice, 
single-spin-flip dynamics will only explore an exponen- 
tially small fraction of the possible configurations. 

So we are forced to turn to a cluster update algorithm 
to simulate this model. The algorithm we use in this pa- 
per is the zero-temperature limit of the Wang-Swendsen- 
Kotecky (WSK) cluster algorithm |l|,|l|l for Potts anti- 
ferromagnets, which is defined as follows. At each step in 
the simulation we choose two of the four colors on the lat- 
tice, identify all connected clusters containing only these 
two, and, in each cluster independently, either switch the 
two colors or leave them untouched with probability ^ 
Clearly this preserves the property that neighbors have 
differing states, and, as we will show, it is computation- 
ally quite efficient. 

At finite temperature the WSK algorithm is trivially 
ergodic, as demonstrated by Wang et al. in their original 
paper. With a little more effort it is also possible to 
prove that it is ergodic at T = for all g > 2 on the 
square lattice |l^] or more generally on bipartite lattices 
|O,E0[. The triangular lattice however is not bipartite, 
and in fact the algorithm is known not to be ergodic 
at T = on triangular lattices with toroidal boundary 
conditions, at least for certain lattice dimensions iQ. In 
this section we make use of our height representation to 
prove that the algorithm is ergodic for certain other types 
of boundary conditions. 

Previous proofs of the ergodicity of the WSK algo- 
rithm at zero temperature have relied on defining a spe- 
cific "target configuration" of the spins on the lattice 
and demonstrating that any given configuration can be 
transformed into this one in a finite number of reversible 
moves. For bipartite lattices the target configurations 
used have been checkerboard colorings. Here we use the 
same approach to prove the algorithm ergodic on the tri- 
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FIG. 10. A three-coloring of the triangular lattice. There 
are six such colorings for each choice of three colors. 



angular lattice but, since the lattice is not bipartite and 
therefore does not permit checkerboard colorings, we use 
instead a three-coloring of the lattice as our target con- 
figuration. 

There are six possible three-colorings of the triangular 
lattice for each choice of 3 of the 4 colors. We illustrate 
one of them in Figure |l^. We define a domain to be a 
connected set of sites whose colors coincide with one of 
these three-colorings. The three colors in such a domain 
fall on three triangular sublattices with lattice parameter 
-s/S times that of the fundamental lattice. Our goal is to 
add new sites to the domain, singly or in groups, until 
the domain fills the entire lattice. If we can do this with 
reversible cluster moves from any given starting state, 
then our algorithm must be ergodic, since we can get 
from any state to the target and then from the target 
to any other state. (It makes no difference which three- 
coloring we take as our target configuration since we can 
get from any one to any other in at most three steps of the 
WSK algorithm — one to get the three colors right and at 
most two to put the colors on the correct sublattices.) 

Our approach is as follows. We choose a site which lies 
outside our domain, but is adjacent to it. The color of 
this site, call it a, differs from the color b of the sites inside 
the domain on the same sublattice. If we can switch 
the colors a and b everywhere inside the domain, while 
leaving the color of the new site unchanged, it will now 
match the other sites on that sublattice and we will have 
added it to the domain. (Note that the domain now has a 
new three-coloring, resulting from our switching a and b.) 
This can be done trivially in a single Monte Carlo move 
if we make the right choice of two colors for the move, 
but there is a catch. The problem is that there may be 
some cluster connecting the new site to a (possibly quite 
distant) site in the domain via sites of colors a and b. If 
such a cluster exists, then the new site will get flipped 
whenever we change colors within the domain, so that its 
color will always differ from that of the sites on the same 
sublattice in the domain. An example of such a situation 
is shown in Figure O. 

Happily, we can invoke the height representation to 
prove that any state in which such a path exists must 
contain at least one defect, and hence that such paths 
cannot exist in a ground state of the system. To see this 
we consider a closed loop of sites formed by the path out- 
side the domain, completed with any path of our choice 
within the domain. To make things concrete, let us sup- 
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FIG. 11. A configuration which defeats the WSK algo- 
rithm. We would like to grow the domain indicated by the 
dotted line by adding to it the spin with state 4 immediately 
above it and propose to do this by changing each of the 3s in 
the domain into a 4. But a path of 3s and 4s connects the 
domain to the new site, so it will be changed to a 3 when we 
do this. In the text, we show that configurations of this kind 
are forbidden within the ground-state ensemble because any 
contractible loop of this kind must contain a defect. 



pose that a and b are the states 4 and 3 as in Figure |T^. 
The portion of our closed path outside the domain con- 
sists by definition of sites of only these two colors and 
hence the change in height Ah from one end to the other 
must be a multiple of, in this case, a. Within the domain, 
on the other hand, the heights on all sites belonging to 
a particular sublattice are the same, so that Ah for the 
portion inside the domain is just equal to the change in 
height resulting when we change a 3 in the domain to 
a 4, or /3 — 7 in this case, which is perpendicular to a. 
In general. Ah for changing some site in a three-coloring 
from one color to another is perpendicular to the Ah be- 
tween neighboring sites of those colors. This means that 
the sum of the Ahs for the two parts of the loop cannot 
be zero, and hence the Burgers vector around the loop is 
non-zero and the loop must contain a defect. 

Unfortunately, this does not quite prove the ergodicity 
of the algorithm. Certainly a configuration like that in 
Figure ^can be ruled out, because there must be at least 
one defect within the loop, and hence the lattice cannot 
be in a ground state. (The reader might like to popu- 
late the interior of the loop with spins just to check that 
this is indeed true.) However, there is another possibil- 
ity. If we have some form of periodic boundary condition 
on our lattice, it may be possible for the loop to go off 
one side of the lattice and come back on another and re- 
join the domain that way. It turns out that such a loop 
can have a non-zero Burgers vector even when the lattice 
is in a ground state. In essence, the lattice possesses a 
non-localized defect, without there being a defect any- 
where in particular. The crucial difference between the 
two types of loops is that in the first case the loop is 
contractible, meaning that it can be shrunk to a point 
by shifting it one plaquet at a time. For a loop of this 
type, the situation depicted in Figure |ll| applies, and the 
arguments given above are correct. If however the loop is 
non-contractible, then it is possible for there to be a non- 
localized defect and we cannot prove the ergodicity of the 
algorithm. To give an example, we show in Figure n3 
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FIG. 12. A configuration on the torus with nonzero Burg- 
ers vector around its non-contractible loops. The domain de- 
noted by the line cannot grow monotonically under the WSK 
algorithm. 



a configuration of the model on a lattice with toroidal 
boundary conditions. This lattice has two fundamental 
non-contractible loops, one wrapping around the bound- 
ary conditions horizontally, and one vertically. For the 
configuration shown, these two loops have Burgers vec- 
tors of 2f3 + 47 and 4/3 -I- 2j respectively, even though 
there are no localized defects. In fact, the domain inside 
the dotted line shows how the algorithm can fail. Any 
attempt to add a new site to the domain by exchanging 
two colors inside the domain will change the color of the 
new site as well, because the new site is connected to 
the domain by paths which wrap around the boundary 
conditions. This does not actually prove that the WSK 
algorithm is not ergodic on this lattice, only that we can- 
not prove it to be so using arguments of the type given 
here in which domains gain sites but do not lose them. 
However, Salas and Sokal iQ report that they found a 
configuration on a 6 x 6 torus which can only be trans- 
formed into a small number of others, thus showing that 
the algorithm is not in general ergodic on the torus. Huse 
and Rutenberg, in their study of the three-state model on 
the Kagome lattice, noted a similar lack of ergodicity in a 
loop-flipping algorithm on the torus |1^. It seems, there- 
fore, that it would be imprudent to conduct simulations 
solely on a toroidal or similar lattice. 

On lattices with no non-contractible loops, however, 
the algorithm is, by the arguments given above, ergodic. 
Examples of such lattices include the infinite lattice and 
any finite lattice with free boundary conditions, where 
opposite edges of the lattice have no connection to each 
other. Unfortunately, the first of these is impractical for 
computer simulations, and the second suffers rather dra- 
matically from finite-size effects. A better solution is to 
perform our simulation on a lattice with periodic bound- 
ary conditions, but with a topology chosen so that there 
are no non-contractible loops. The simplest example is 
the sphere. Because its Euler characteristic is 2, there is 
no way to cover a sphere with a regular triangular lattice; 
some sites have to have less than six nearest neighbors. 
Fortunately, our proof works not just for the triangular 
lattice, but for any lattice where both the height repre- 
sentation and the target three-coloring can be defined. 
As we showed in the previous section, the height repre- 
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FIG. 13. Left: an octahedral lattice with L = 4. The WSK 
algorithm is ergodic on this lattice. Center: a three-coloring 
of this lattice around one of the vertices of the octahedron 
and the definition of the orientation vector t along the bonds. 
Right: the structure is that of a triangular lattice everywhere 
except at the vertices of the octahedron. 
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FIG. 14. An L = 3 lattice whose topology is that of a 
projective plane. Pairs of sites connected by dotted lines are 
identified. Note the consistency of the three-coloring. 



sentation can be defined for any planar lattice with tri- 
angular plaquets and even coordination number, so that 
we can place vectors t on the bonds oriented so that they 
go only clockwise or only counter-clockwise around each 
plaquet. We can also define a three-coloring on such a 
lattice by having colors cycle upward or downward along 
t. Thus, we can for instance perform a simulation on a 
lattice tiling the surface of an octahedron where 6 sites 
have four nearest neighbors, and the WSK algorithm will 
be ergodic. In Figure |l^ we depict one such lattice and 
illustrate how the three-coloring behaves near its corners. 
Another possible topology is the projective plane — a 
hemisphere with diametrically opposing points on the 
equator identified. To create a lattice with this topol- 
ogy, we simply take the upper half of the octahedron 
depicted in Figure O and identify sites along its square 
equator as shown in Figure p^. Since the projective plane 
has only half as much curvature as the sphere (its Euler 
characteristic is 1) there are now just three sites with four 



neighbors rather than six. Thus this lattice is propor- 
tionately flatter than the sphere with the same number 
of sites. Note that the three-coloring is well-defined for 
any lattice size. 

Although the projective plane does possess non- 
contractible loops, every such loop has the property that 
going around it twice makes it contractible. Formally, 
the fundamental group Hi for the space is Z2, the inte- 
gers mod 2. This means that the Burgers vector of any 
such loop must satisfy A/i -I- A/i = 0, so Ah = 0. (Kawa- 
mura and Miyashita [p2| show that the Heisenberg an- 
tiferromagnet on the triangular lattice has defects with 
a Z2 charge, so there are spin systems in which Ah can 
be nonzero even though 2Ah — 0. This however is not 
the case with the present system.) Thus non-contractible 
loops are tolerable as long as they have finite order, i.e., as 
long as their n**^ multiple is contractible for some finite 
n. The sphere and projective plane are the only finite 
two-dimensional manifolds satisfying this condition, since 
both tori and Klein bottles have non-contractible loops 
of infinite order. 

While Monte Carlo calculations performed on spheres 
and projective planes may seem outlandish, they are im- 
portant in the present case in order to rule out the possi- 
bility that lack of ergodicity is introducing a bias into our 
results. The drawback is that we are forced to introduce 
a small number of atypical sites into the lattice (those 
with only four neighbors) which, for example, makes cal- 
culation of spatial correlation functions more difficult. In 
this paper we strike a compromise by performing some 
(probably non-ergodic) simulations on toroidal lattices, 
which give excellent statistics for correlation functions 
and other quantities, and some on the projective plane, 
which give poorer statistics, but are more trustworthy. In 
fact, we find that there are no physical measurements for 
which the two topologies disagree, so it is possible that 
the simulations on the torus are "sufficiently ergodic" for 
our purposes, although we cannot guarantee this. 

Just as with the height representation, our proof of 
ergodicity will work for other lattices as well, whenever 
we can give each bond a direction such that the bonds 
around each triangular plaquet go either only clockwise 
or only counter-clockwise. It can also be used for some 
lattices which have a mixture of ferromagnetic and anti- 
ferromagnetic bonds, but which have no frustration. At 
T = vertices with a ferromagnetic bond between them 
can simply be identified, since they must have the same 
state. Then the WSK algorithm is ergodic on such lat- 
tices if the resulting graph fits the conditions above. 



VII. RESULTS OF MONTE CARLO 
SIMULATIONS 

We have performed extensive simulations of the four- 
state triangular model using the WSK algorithm on lat- 
tices having the topology both of the projective plane (for 
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FIG. 15. Correlation times for the WSK algorithm on the 
projective plane (circles) and the torus (squares) measured in 
Monte Carlo steps. Both lattices give a dynamic exponent of 
z = 0.74 ±0.02. 



which the algorithm is definitely ergodic) and of the torus 
(for which it probably is not). On the projective plane 
we simulated systems with linear dimension L equal to 
a power of two from L = 4 up to L = 1024. On the 
torus the three-coloring of the lattice is only well defined 
for L a multiple of three, so we simulated systems with 
L = 6, 12, 24 . . . 768. In each case simulations ran up to 
one million Monte Carlo steps for the largest lattices. In 
order to allow accurate error estimation, and also to ex- 
amine the efficiency of the algorithm, wc first measured 
the (exponential) correlation time r in Monte Carlo steps 
as a function of system size. Figure Ha shows the results 
for both topologies. In each case the dynamic exponent 
z, defined by r ~ L^, takes the value 0.74 ± 0.02, which 
is comparable with values for cluster algorithms for fer- 
romagnetic Potts models p^. Correlation times for the 
projective plane are about a factor of 1.3 larger than 
those for the torus with the same value of L, which is 
presumably because the projective plane has more sites 
on it. (There are L^ sites on a torus of linear dimen- 
sion L, but 2L^ -f 1 on the projective plane.) 

In practical terms, since a single step of the WSK al- 
gorithm updates a number of spins which scales with the 
area L'' of the lattice, the CPU time taken to generate a 
given number of independent lattice configurations scales 
as L''+^ ~ L"/4. 

In order to measure the critical exponents for the 
model, we chose two different definitions for the order 
parameter and measured two-point correlations and fluc- 
tuations for each of them. From these results we ex- 
tracted values for the correlation exponent 77 by direct 
measurement, and for the susceptibility exponent ^/i^ by 
finite size scaling. These exponents are not independent; 
we expect them to be related by the Fisher scaling law 
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FIG. 16. The susceptibility for the staggered magnetiza- 
tion on the projective plane as a function of lattice size L. 
The dotted line is a fit to the last five points and gives 
■y/u = 0.71 ± 0.02. Theory predicts j/iy = ^■ 



0.85 



0.80 - 



■^ 0.75 

-a 



S 0.70 



0.65 



0.60 



10 



100 



1000 



FIG. 17. Estimates of 7/1/ obtained from fits to simulation 
results for the staggered magnetic susceptibility on lattices of 
size ranging from L up to 1024, as a function of L. These 
tend downward with increasing L, presumably as a result of 
logarithmic corrections to the scaling form. 



7/^^-2-77 ||,| 
defined as follows: 



15|. Our two order parameters were 



1. The simplest magnetization measure is just rrik = 
Nk — jN where N^ is the number of spins on the 
lattice in spin state k. The two-point connected 
correlation for this magnetization, averaged over fc, 
is then Gc{i,j) — Sg^sj — j and the susceptibility is 

2. We also examined the "staggered magnetization" 
defined by Husc and Rutenberg |1^ for the three- 
coloring model on the hexagonal lattice, to allow 
direct comparison with their results. In the con- 
text of the present model, this magnetization may 
be thought of as a complex number ^^ ciri where 
the sum is over bonds on the triangular lattice, ci 
represents the color a, /3 or 7 of bond I as de- 
fined in Equation (M), and n is a static reference 
coloring corresponding to the color of bond I in 
any of the 24 possible three-colorings of the lat- 
tice. The staggered magnetization mi on a site i 
can be defined as the sum over the bonds / con- 
nected to that site and the two-point correlation is 
then given by Gc{i,j) — rriim* and the susceptibil- 
ity by x = X],m,m*. 

The staggered magnetization is in fact slightly the eas- 
ier of these two to analyze, so we examine this case first. 
Huse and Rutenberg [^ pointed out that the staggered 
magnetization has a wavelength v3 on the height lat- 
tice. Using Equation (Hq) to transform from Cartesian 
coordinates to triangular ones. Equation (HS) then gives 
77 = 3 and j/v = |. Figure |l^ shows our simulation 
results for the susceptibility on the projective plane. A 



least-squares fit gives ^ /v = 0.71 ± 0.02, which is a lit- 
tle greater than the expected value, but, as Figure H^ 
shows, there is a clear downward trend in the value with 
increasing system size. Huse and Rutenberg saw simi- 
lar behavior in their bond-coloring model (in fact their 
value for ^ /v is almost exactly the same as ours) and 
they attributed it to logarithmic corrections to the scal- 
ing forms which arise because the height representation 
is at its roughening transition (see Section IV). 

We have also measured spatial correlations in the stag- 
gered magnetization on both the torus and the projective 
plane (Figures |l^ and |^, respectively). The results on 
the projective plane have larger statistical errors than 
those on the torus, because of the need to stay well away 
from sites with local curvature in performing the calcu- 
lations. Fits to the data yield values of 77 = 1.27 ± 0.01 
on the torus and 77 = 1.28 ± 0.09 on the projective plane, 
in reasonable agreement with theoretical predictions. 

The situation is a little more complicated for our other 
definition of magnetization. Looking at Figure M, we see 
that the mapping from heights to spin states has wave- 
length 2 on the height lattice, which implies that both the 
correlation exponent for the spins and the corresponding 
susceptibility exponent should be equal to 1. Figure |20| 
shows our data for the susceptibility on the projective 
plane. A least-squares fit gives ^/v = 0.84 ± 0.01, which 
is some way from the theoretical prediction, but, as Fig- 
ure pi] shows, the value is increasing steadily with system 
size, so the discrepancy is again probably due to logarith- 
mic corrections. 

However, when we look at the spin-spin correlation 
function. Figure g2[ we find that 77 = 0.35 ± 0.01, which 
is nowhere near 1. Our values for the exponents thus 
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FIG. 18. The connected correlation Gc(f) for tiie staggered 
magnetization on the torus. The best power-law fit is indi- 
cated by the dotted hne and gives r] = 1.27 ± 0.01. Theory 
predicts 7)= \. 
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FIG. 20. The susceptibility for the (un-staggered) magneti- 
zation on the projective plane for systems of linear dimension 
L up to 1024. The dotted line is a fit to the last five points 
and gives -^jv = 0.84 ± 0.01, while theory predicts 7/1^ = 1. 
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FIG. 19. The connected correlation Gc{r) for the staggered 
magnetization on the projective plane. While the data are not 
as good as those for the torus in Figure PlSl we get a compat- 
ible value of 7; = 1.28 ± 0.09 for the correlation exponent. 



FIG. 21. Estimates of 7/1/ obtained from fits to simulation 
results for the un-staggered magnetic susceptibility on lattices 
of size ranging from L up to 1024, as a function of L. These 
tend upward with increasing L. 
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FIG. 22. The connected correlation Gc{r) for the 
(un-staggered) magnetization on the torus. The best 
power-law fit is indicated by the dotted line and gives 
T] = 0.35 ± 0.01. Theory predicts n ^ k- 
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FIG. 23. The sublattice of possible height differences be- 
tween two sites along a principal lattice direction, both with 
spins in state 1, whose separation on the lattice is a multiple 
of 3 times the lattice parameter. The dashed path 2/3 -I- 4q: 
is allowed, but the dotted path 27 + 2a — 2/3 is not, since 
it contains terms of different signs. The wavelength of the 
sublattice is 2\/3, which implies that J? = | for the spin-spin 
correlation function. 



appear to violate the Fisher law. The explanation of this 
result is as follows. 

The correlation function shown in Figure E2 is mea- 
sured, naturally enough, along one of the three principal 
directions on the toroidal lattice. If two sites lie along 
such a direction and their spins are in the same state, 
then their heights are constrained to particular sublat- 
tices of the height lattice. For instance, suppose that the 
distance between them is a multiple of 3, and that one 
site has state 1 and height (0, 0). Then the height of the 
other site is a sum of multiples of three times a, (3 and 
7, all with the same sign since t is constant along any 
of the lattice directions. Since q; + /3 + 7 = 0, we can 
cancel terms in threes until only two kinds of terms are 
left, say as and (3s. Since the number of these is still a 
multiple of three, and since a and [3 (or more precisely, 
the unit vectors in those directions on the height lattice) 
have a component of ^ along the hi axis, the only heights 
we can end up with are ones with hi = 3fc/2 for some 
integer k. The set of such heights which correspond to 
spin state 1 is shown in Figure ^. More generally, once 
we choose the distance between two sites and the state of 
the spin on one of them, the only heights the other site 
can have and still be in the same spin state are those on 
one of the three sublattices corresponding to that state. 
The wavelength of each such sublattice is 2\/3, and Equa- 
tion na then gives 77 = 1/3, in good agreement with our 
simulation results. 

This phenomenon, in which some constraint gives an 
operator a wavelength on the height lattice longer than 
that for the corresponding susceptibility, may apply to 
other operators and systems like this one. Therefore, it 
appears that the Fisher scaling relation cannot always be 



applied in a straightforward fashion. 



VIII. CONCLUSIONS 

In this paper we have studied the four-state Potts an- 
tiferromagnet on the triangular lattice at zero tempera- 
ture. By examining the spectrum of defect types which 
can appear in the model and identifying the conserva- 
tion laws which govern their interactions, we have been 
able to define a Burgers vector for the model and thus 
show that the ground-state ensemble has a well-defined 
height representation. The height is, in this case, two- 
dimensional and may be the simplest example to date of 
a vector height. 

Using the height representation, we have been able to 
demonstrate a number of results. We have shown that 
the model is equivalent to a three-state Potts antiferro- 
magnet on the bonds of either the triangular or hexagonal 
lattice, or on the sites of the Kagome lattice, that pairs 
of defects feel entropic forces between them in proportion 
to the dot product of their topological charges, and that 
the spin-spin correlations in the ground-state ensemble 
must decay algebraically at large distances. 

We have calculated exact values for a variety of crit- 
ical exponents. The scaling exponent rj for the spin- 
spin correlation function is of particular interest because 
the wavelength of the spin operator on the height lattice 
turns out to be longer than the fundamental periodic- 
ity of the height-to-spin mapping due to a constraint on 
paths connecting sites along a principal lattice direction. 
This gives a value of 77 = i even though the correspond- 
ing susceptibility exponent 7/r/ = 1 — an apparent viola- 
tion of the Fisher scaling relation. 
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We have also used the height representation to prove 
for the first time that the Wang-Swendsen-Kotecky clus- 
ter Monte Carlo algorithm is ergodic for the four-state 
model on the triangular lattice. Our proof however re- 
quires that the lattice have a topology which possesses no 
non-contractiblc loops of infinite order, and this means 
that simulations on the torus (which has such loops) are 
probably not ergodic, calling previous simulations of this 
and related models into question. Simulations on lat- 
tices with free boundary conditions or lattices with the 
topology of the sphere or the projective plane are, on the 
other hand, provably ergodic, and we have performed 
extensive simulations on the projective plane. We find 
reasonable agreement between the values of the critical 
exponents measured in these simulations and both theo- 
retical values and values from previous numerical studies. 
There are, however, significant logarithmic corrections to 
scaling associated with the fact that the height model is 
exactly at its roughening transition, and this means that 
we have to go to extremely large lattice sizes to see the 
expected behavior. 
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Note added. During the course of these investiga- 
tions, we learned of an unpublished manuscript by Hen- 
ley which also shows that the four-state antiferromagnet 
studied here is equivalent to the three-state one on the 
Kagome lattice. 
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